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* SUMMARY 

The practicability of using a stretched cable for a rocket 
launching device is dependent primarily upon the transverse move- 
ment of the rocket and the stresses involved in the system. This 
requires analysis of the effects of the moving rocket mass on the 
cable from the instant of contact, and the problem reduces to onc 
of wave propagation. 

The analysis here is restricted essentially to developing a 
procedure by which the initial rnotion of the mass can be calculated, 
and an expression is obtained which permits determination of the 
deflections. The problem io approached by first assuming that the 
mass exerts a constant transverse force on the cable. Admittedly, 
this is a simplifying assumption, and the result is not valid for the 
instant the mass hits the cable or for a short time thereafter be- 
cause the inertia of the mass is not considered. 

Next, the problem is solved by taking into account the dy- 
namics of the mass, and the solution reveals that the path of the 
actual mass deviates from the force path by as much as twenty per- 
cent during the initial motion but svon returns to the force path. 
The mass does not exert its full force on the cable at the instant of 
contact, but comes down on the cable with full force a short time 
later. Since maximum stresses oh the system occur at this time, 
this factor is an important result of the analysis. Further, itis 
shown that increased velocity of the mass increases the deviation 
from the force path, and increased mass lengthens the time of 


return to the force path. 
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ZI 
I. BENIRODUCTION 

The problem investigated here results from an original 
inquiry by the U. S. Navy as to the feasibility of using a stretched 
cable for a rocket launching device. Recent interest in the subjcct 
has extended greatly the range of os gible applications, and since 
the analysis of the physical phenomena involved in a mass moving 
on a stretched cable is general in nature, it need not be restricted 
to the launching of rockets. 

С. ЕК. DerPrima and Е. ЕБ. e have investigated thc 
problem of a force moving on a stretched cable. Their results in 
the supercritical case (i.e. where the force travels faster than 
the — of propagation of a disturbance on the cable) show that 
although the cable distorts under the influence oí the force, the 
path of the force is unaffected and all distortions occur behind the 
moving force. In the subcritical case, the path of the force is af- 
fected but the initial conditione arc unrealistic because a force 
and not a mass was considered. Of interest then is the subcriti- 
cal case of a mass traveling with conatant velocity across a 
stretched cable. By considering the subcritical case, the analy- 
sis will hold for slower and heavier masses, and the range of 
application can be extended. 

For the purposes of this work, a portion of this phase of 
the problem is investigated: i.e. the initial motion of a mass 
traveling with constant subcritical velocity on a stretched cable. 
The procedure will be to derive the partial differential equation 


which describes the physical phenomena, solve the equation, and 















* wae 
П. DEVELOPMENT OF THE PROBLEM 


For the system, 
assume a stretched string 


between two fixed points 





and consider a mass mov- 
ing across it with constant 
Fig. 1 velocity v (Fig. 1). 
First, consider the partial differential equation for the trans- 
verse motion of a string 
=~ distorted into some curve 
u(x,t), (Fig. 2). Applica- 
u tion of Newton's second 
Fig. 2 law leads to the wave equa- 


(2) 


tion in one dimension 


es zT Pau (1) 
* NET 
8t ox 
where p ig mass per unit length of string, T is tension in string, 
2 
and u(x,t) is the displacement of the string. The term pos is 
| 2 дї 
да 


the inertia force due to string mass and Ty is force due to ten- 
Ox 
sion in the string. 
If an arbitrary force F(x,t) is applied to the string, Equa- 


tion (1) becomes 


Pa". TT 
Ox c Ot 
where c* = 5. 








+ -4- 
If the force is exerted by a mass m on the string at a def- 
inite point, x distance from the origin (Fig. 3), and the mass has 
a displacement W{t), then 
the force F applied at the 
point of the mass may be 


calculated by applying 


Newton's laws to the masa. 





or 
F = -їп (—- - 2) (3) 


where mg is the force due to gravity. 

The acceleration of the mass may be expressed in terms of 
u, the displacement of the string. Consider the mass moving across 
the string with uniform velocity v. Taking an instantaneous look at 
a particular point P, the vertical velocity of the mass is = 
However, the verti- 
cal velocity due to the 
velocity v along the 
string must also be 
considered. Looking 
at Fig. 4, (greatly 
exaggerated because 
9 actually is very 


Fig. 4 small), 





e 


. ¿5 
Vertical component of v = v sin Ө 


and replacing sin 6 by tan 6, 


A 
Vertical component of V = v tan © = v= 
Au 8u 
Passing to the limit, lim — c 
Ax>0 Ах Өх 


and therefore 
Vertical component oí v = #7 
dx 


Thus, the total vertical velocity in terrns of u(x,t) is the 


du д 
sum of two terms: FẸ tv 3 It was agreed before the vere 


tical displacement of the mass at a fixed point is Wit). Therefore, 


if u is fixed at a point determined by x,t, 


aw du 
wel, 


d ,dW 
At Set oe Get’ oe ae 











ey ot, 2 au ‚au, ‚du 
~ 8xót 2 2 59% 
Ox Ot 
or 
ағу bA 95^u 2 844 
4777 EE У 
dt ot Ox 
x,t 
and introducing this into Equation (3) 
2 2 2 
F = -m (2 + 20 put - с) (4) 
ot ox 
x,t Ket x,t 





Е -6- 
Since at t = 0, u = 0, and the mass moves with a constant 
velocity v, the position x of the mass is fixed by x = vt. Using this 
notation to denote a particular point and combining Equations (2) 


and (4) gives 





u 1 atu m 8u atu 2 9^u 
Ox c t Ət 9х x=vt 


For all values of x other than x = vt, Equation (5) gives the 


homogeneous wave equation, i.e. 


for x # vt 


and therefore it is convenient to introduce the Dirac delta func- 


tion! 3) defined as follows: 


vt+4 € 
O for x e vt 
d (x-vt) z (x-vt)dx = 1 
СО for x = vt 
vt-c 


or saying the same thing 


t to 
[ites wane (6) 


Introducing this function, Equation (5) becomes 


2 2 2 д 2 
— E es +, 2 E qe x Fr + ye 25 - g) S (x-vt) (7) 
x С t t x 


This is a linear inhomogeneous partial differential equation of the 
hyperbolic type which is complicated by the fact that the coefficiente 
are variables. The problem at hand is to solve this equation, and a 


dimensionless forn. will be convenient. 





p 


Define the dimensionless parameters 


u 

UE == . ing 
і-Т уЕ д 
ct V с“ 
Ды: ARE = тт 


where L is the length of tlıe string between the two fixed points. 
First, note ДЖ (х-уї)ах = 1 and since dx has the dimension 
of length, this irnplies the delta function has the dimension u 
This indicates L ۳ (x-vt)] is dimensionless, and dimensionless 
notation for the delta function is Y (5-47) *. Therefore 
L [S (x-vt)] corresponds to а (Є AT). 4 
Multiplying the left side by М ‚the right side by с, ; 


and noting T z pc^, Equation (7) can be written 


e^) ez) : a) or) ” a(t) 
aE)? aS)? Te Fe NU Es 
D 
+ (У) P eu L S(x-vt) 


AF r) 


от 


ay Oy oa fi 8 OY да we Moh sis. 2 008 
os? ar? ¬ A ar? sör os? A (8) 
It should be noted that the initial conditions imposed by the fixed 
end point of the stretched string are v = A =Oat 7 =0. 

Two solutions will be considered in the analysis of the prob- 
lem: first, the solution by assuming a constant transverse force on 
the string; and second, the solution of Equation (8) which takes into 


account the effects of a mass on the string. The constant transverse 


* . xwv ct Ет » 
Se = vt can be written == — — which becomes <` AT. 
L c L $ 7 


o 





-8 ~ 


force equation is obtained by making 8 = 0 in Equation (8) and is 


8*v 9%, 


-= -AÀA /(є-дТ) 
og? Я SA (9) 


92 
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ПІ, SOLUTION OF THZ PROBLEM WITH 


CONSTANT TRANSVERSE FORCE 


Courant ®) shows that the solution of Equation (9), with the 
given initial conditions, can be written 
Т SECO) 
v(£,t)s 4 Sf (x-A t) dxdt (10) 
tz0 xzf-(t-t) 
where g,” arc the dimensionless coordinates of the physical sys- 
tern, and «x,t are used as variables of integration. It should be em-e 
phasized that the variables of integration x,t as used in Equation (10) 
are not the same as the x,t coordinates of the physical systern used 
in the original development of the problem. 

Since interest is focused on evaluating the solution on the path 
of the moving force defined by $ =ar , у(6,7) can be written 
v{iur,r)= WA ^) where W * ) represents the displacement of the 
force. Equation (10) becomes 

Пе) 
Wor) =4 С(х-д 9) dxdt (11) 
t=0 x=¢ -(7-t) 

The evaluation of the integrals of Equation (11) can be accom- 
plished readily by geometric interpretation and the use of character- 
istics. Since the slope of the characteristics in the dimensionless 

§-?t plane are + l and the subcritical case (v < c) is to be consid- 
ered, the diagram (Fig. 5) can be drawn. The path of the force, 


labeled I, is given by E ZAT and is drawn with a slope 14 5, 40. 





-10- 





Гір. 5 


The initial condition imposed by the physical fact that the 
string is fixed at x = 0 may be taken into account as follows. If 
the string is assumed infinite in length, then the efíect of this fixed 
end is simulated by assuming an equal force starting at x = 0 and 
traveling in the negative direction with velocity v. Similarly, the 
boundary condition imposed by the string being fixed at x = L can 
be simulated by an equal force starting at x = 2L and traveling in 
the negative direction with velocity v. The paths of these fictitious 
forces are labeled -F. As will be shown, these patlıs have a neg- 
ative effect in the geometric process for evaluating the integrals. 
Further, it should be noted that the distance 1. in the dimension- 
less plane is unity. 

The arbitrary points (€ , 7" ) are picked on the path of tue 


force and triangles of integration (commonly referred to as 





-11- 
domains cf dependence) are defined by drawing characteristics 
through them. From Fig. 5, it immediately is evident that two 
cases must be considered if a solution is to be obtained along the 
entire length of the string. At the point СЫ 7”) „ the force en- 
counters the wave which has been reflected from the fixed end, 
and the solution at point (¢/%) , good for O<T<t”, will be 
different than the solution at point (( 5 Т") , good for т" < Т 2 
For the purposes of this work, only the initial motion of the force 
is considered and therefore detailed analysis is presented for just 
the first case, i.e. where 0 44 < T*. Proceeding now to the eval- 
uation of Equation (11), 


it will be convenient to 


refer to the simplified 





diagram, Fig. 6. 
(51,7%) 


O 
Fig. 6- 


From the definition of the delta function, 


g+ (7 ~-t) O for x #ut 
JIx-ut) dx = 1 for x = ut 
х-6-(7-% -lforx=-yzt 


and since the value of the solution is dependent on the given values 
only іп the interior of the triangle of integration, this integral yields 
a value of + lat, and only at, each point on the line $=47 , and 
a value of - 1 at, and only at, each point on the line 6--4” . Tnis 
may be thought of as a series of positive unit impulses along є = 47 


and negative unit impulses along $=-4P  , the efíects cancelling 








e -12- 
in the region 0 to 7 . Now taking into consideration the summing 
Ж 
process indicated by d. dt, the double integral of Equation (11) yields 
¿(20 


nothing more than T- A. It can be shown that 


t-t = ADIs = AP 2 


1+ а lta 
Therefore 
_ Ад 
W(t) = 9л (12) 
It is evident that a 
м (7) ' 
plot of displacement versus 
Slope = A, 

\ time would give a straight 

И line, arbitrarily drawn in 

Fig. 7 Fig. 7. It should be noted 


that positive displacements are directed down or earthward in the 
physical sense. 
Investigation of the physical results of this solution reveals 
that if a force F = © is moving along the solid foundation with con- 
stant velocity v at t < 0 and 
hits the string at t = 0, the 
| string deflects instantaneously, 
| and the force takes on a cone 
stant vertical velocity (Fig. 8). 
Therefore, the vertical veloc- 
ity jumps from zero to some 
finite value and the accelera- 


Fig. 2 tion is infinite att = 0. This 


^ 


a -15- 
of course is physically impossible. 

The reason for such a result is obvious. Equation (11) gives 
the solution for a constant transverse force moving across the string 
and not a mass which is physically the case. That is, Equation (11) 
ignores the dynamics of a mass. Intuitively, it is evident that a mass 
moving along the solid foundation and hitting the string would give a 

path probably something 
like that of Fig. 9. The 
inertia of the mass would 
prevent in effect the infin- 
ite vertical acceleration. 
The solution of Equation 


Fig. 9 (8) which takes into account 


the effects of a mass should correct this discrepancy. 





0 -14- 
IV. THE DIFFERENTIAL DIFFERENCE EQUATION 
FOR THE MOTION OF THE MASS 


Rewrite the complete physical equation (i.e. Equation (8) ) 


ay ay. 1-80 a al OO Sigan (13) 
052 ap? ar 9597 ee 


The solution in integral form where x,t are the variables 


of integration is 


¿+7 -t) 


v($,T)- a ç E ^ b^ N ov ) {(x-pt)dxdt (14) 
е 9хәӘї Ox 
t=0 x=€- Pu. 





= a | v 
A 

\ 
D 


N 


N ыы 
\ Y ON 


Fig. 10 
It will be remembered that the integral with respect to x 
is zero except at each point on the path of the mass where it yields 
unity. Further, the values of this integral cancel in the region 


0 to Y (Fig. 10). It can be shown easily that 








| 


le 
= i) r (15) 
and since the path of the mass is defined by £ =w Y , Equation (12) 


can be written 


2 
A UN rap EA dt (15) 


‘Using the i А W(r) = v( mr, Y ), it io clear that Equa- 


tion (16) becomes 
YP 


ab gia ug. M dt (17) 
t = at 


where the system parameter a= 4 


Integrating Equation (17) gives 


W(e) 2 30-2 | eget pen) for 0« t « t" (18) 


This is the differential difference equation which can be used to con- 
nect the solutions at two values of time, ar and f. It should be noted 
that the validity of this equation ie limited to 0 « T « t^ 

In a similar rnanner, a differential difference equation can be 


developed for that portion of the string where a <д , (Fig. 11). 


N 
(ETT) / j -M oS, N 


I—II اجس‎ 


Fig. li 





alga 
* 
It is simple to show 





lis ы 
ті: (207 (12) 
1+ 2 
T, - e m - TA (20) 
2 
ues dex (21) 


Writing Equation (14) for this casc, it is clear that 


"m T, 
2 2 
t=T, t= f3 


Integrating and using Equations (19), (20) and (21) gives the 


differential difference equation good for MUT « ni 


сез = 2240-77 A - aut. ato, alt ا‎ (22) 


This equation will permit evaluation of the displacements of the mass 


over the last portion of the string. 








Е ~l م‎ 
YV. SOLUTION FOR THE INITIAL. MOTION 


Rewrite /;quatioa (12) 
31574, ау; 
W(t) = E (1-0) ~ E Es + ze un | (23) 


This differential difference equation connects the value of the dis- 
placeınent of the глаза W(T ) at two values of time within the liinits 
0 <r<t”" . It should be noted that a 7 1 and conscquently or 
defines times less than T. Bearing these facts in mind and using 
a stepwise procedure, this unique equation can be used to derive 
an expression for the initial motion of the mass. The number of 
steps taken will determine the accuracy of the result. Three steps 
will be sufficient for the purpose of illustrating the procedure. 

By using Taylor's series expansion and expanding about zero, 
it is possible to obtain an expression for the solution which will in- 


2 clude a point, say a p 







close to zero on the path 
oí the mass (Vic. 12). К 
the next step is taken to 
point ат, the solutions 
at the two points are Cone 


Region of nected by the differential 


Toylor's Expansion 

difference equation, and 

by inserting the result at 
ne ae Е: 

point a T into “£quation 

Fig. 12 23) and inteprating iron 
© b c 


2% toa’ , the solution at point aT is obtained. Repeating the 








10 


. -). е» 


process will yield the solution at t . It is of interest t note that 

if steps of a Ay are taken where nis an integer, each subsequent 
point is located geometrically in the €-T plane by use of the char- 
acteristics as showa in Fig. 12. Proof of this fact is trivial. The 
actual mechanics of the procedure arecarried out as follows. 


Differentiating Equation (23) 


2 2 
ау”) „A ТЕШ: a’wir) _ | 
dr 2 1 2 "me 3 dlar) i 


and as T —» O0, W(T) 2 W(oT ) 2 W(o9). Therefore 


W'(o) = $ (1-0) - Б (1-а) W (о) 
w оу 7 V (0) + 5 (1-а) 


5 (1-a) 


Using these in the Taylor's series 


W'(o) * à E T 


l-a 


W(T)  W'lo)T f 


and inserting a roore convenient notation waere 
-W'(o) 4 à (1-a) 
A 8 W'(o) Bz 


2 (24) 


W(T) AT + BT” 4 ----- 


It is evident that the accuracy could be increased by taking more 
terms of the series. 


Taking the first step to ar — , Equation (24) gives 
W(a^r ) E Mer ) + Bla? js (25) 


where war ) is interpreted to mean the displacement at the point 


dr А 








. 2192 


The next step to point of ig taken by use of the differentia] 


difference equation. Inserting at for f in Ziquation (23) 


5 Srt + Wart) =2 (l-a lar) +6 er 


2 
d ү a” 
ar) 


Ly» 
Using Equation (25) to obtain «УДат) 
d(a 7) 


E M * W(a^) = 3 (1-0Ҳат) +5 |2+2B0(0 7) | 


This equation now is integrated over the step under consid- 


eration, i.e, from ат toat , to give the result 


дат 2 2 Ser Zar 
e мне” SES (ar) - 4E. ef 


тал А Gar Lar > B 2 E 
-Ae (ат) «2р е tage +Ba p ef (ar) eb 2, 


Inserting the limits and using Equation (25), the expression for War) 


is obtained. 


Repeating the above process to the point T leads to the final 


result 
E. 2 
W(r) = (ET + Ket -к„) E 


FET “FET HE Jw 
8 `9 10 | 26 
UP 2 (26) 
2 > {a` -a +a-1)7 
la KoT | +aK T - Koef 


Ha -1)+ 
* (a, 4 Ку) еЁ tK? + 5 Ky, 








Е -20- 


where 


к 0-06 ве ЕЕ 


7 (a -a+ 1) 
Ka Z ی‎ E в 
P OO RE a1) т 
4 = Е 
K, = Ba K, = a(A-K,) 
Em E S R4 
E, s $(A-X) к. ғ Sê ky 


2 Zu 
Ko = 2Ba Kzll-aK,) + a KEXA-K,) 
27) 
Ко 2 B в к.к, (1-е 5.) - K,ll-a KM KR) ty K3( A. - E) 


2 
kou É (A-K) - B E a 


£,, 24000 ЫП 


By setting e = 0 and renumbering that a = n ‚ Equation 
(26) reduces to 


or 
w(t) = 4-7 (28) 
which is exactly Equation (12), the result of the simple analysis con- 
sidering a constant transveree force. 
The algebraic expression, Equation (26), can be used to solve 
for the displacements of the mass for 0 < T< t* . Selecting the 


fixed parameters to be 


mg = 250 lbs. c = 350 ft/sec. 
1, = 2509 £t. p = 0.1866 slugs /ft. 








a -214 


A and ё can be determined, and calculations were carried out 





for the following representative conditions: (1) A= P W'(o) 2 0; 


адет WHO) m 05 (3) кеў, 0) = FEF ;(4) д=зр Wo) = 0: 
and (5) A= 5 W'(o) = 0 with mass increased tenfold. For conven- 


W(T) 
w. 5 
and (2) are shown in Fig. 13. Condition (5) is shown in Fig. 14, 


ient comparison, plots of - 1 versus Y for conditions (1) 
and the difference in Y scales between the two figures should be 
noted. As should be the case, it was found that conditions (3) and 
(4) were symmetric to conditions (1) and (2), respectively, and 
therefore are not shown. 

Reference to the plots shows that when the dynarnics of 
the mass are considered by solving the complete physical equa- 
tion, the discrepancy of infinite acceleration which appears in the 
result of the constant transverse force analysis is eliminated. In 
fact, the slope oí the displacement curve, W(Y ) versus f , is 
zero at time f = 0. Further, the actual mass path W(T ) returns 
very rapidly to the force path W (r). 

Since 4 = 5 is indicative of a greater horizontal velocity 
than 4 = 5, other parameters being equal, the plots in Fig. 13 re- 
veal that increased velocity gives the actual mass path a larger 
deviation from and slower return to the force path. Also, other 
parameters being equal, the plot in Fig. 14 indicates that an in- 
crease in mass decreases the deviation from the force path and 
increases greatly the time of return to it. 

| Ав was pointed out, greater accuracy can be attained in 
the solution by increasing the number of steps. The result by using 


five steps is given in the appendix. 
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APPENDIX 


p 
5 
= 


A icore accurate 

=pression for the initial 
rnotion of the mass is ob- 
tained by increasing the 
number of differential 
steps taken. The math- 
ematical procedure is the 
same as outlined in the 
text, sirəply being repeated 


more times. If five steps 


> а 
Ww— atr be taken as shown in Fig. 
Pigs. 15 15, the result is 
diu de el)r 
? к ” 
wt) SE E, e + [icr st - 8x, et 
la + -а)7 


3,2 а ar 
+ Bat ta E r AE € кі 523 ef 


pr 
+ | Seat | с 


di -а° 
үк +40 4 E 


a ға 


2 s 
с + Boot “Bon 


Ky Y 4 Ego? “кэз | et 


c 


| а: 


4, 
2 (12-12) 


Gla” =O, 2а „1р 














Г Ne 


2 


| we +a -1) Y 


2 m E. 5 ie 
к лн ну" ү, + Бо) 


KE al} T 


м 


T" a, 37 - (Ko + ap | of 


| Ha -1) D 


5 = (а-а ta -a)T 
+ 5127 + Ko T | e 


Ab -are у 


e 


q rar)? 


^ 
E ла uf t1) ШЕ 
of 


| 
| 
[a spt? t arat ms] 
еее к | 


btts 


= 40 +? M KaT И 


К 
A, +a -l) Y 


25 2 
-fa i kK I кз | 


4 T E c gau +a 2.1) + 
-|e Ki +a Kig “Koo e 


с 
2 fla -a ta -aT 
елү Ks? Su) e 


а” бы el)? 


gla -atta -at 


+ 5 K¿T en as KaT = е 





. -26- 


> N ج‎ 
Жа -atta a") T 


2 qe -ats we -1) 7 
= 2с мәт Mo (С, + Kag) T -4” І эз 4 Ж اوو“‎ | 
> o -a ta -a ta -1)” 
у” > YX 
+ a Esp" + ak, ? Kol 
2 > Ha a ta -ata -1)Y 
+ ж.” IA e 





A +a zo x ~l) T 
ef 


> dete T +a -1)7 
a X ge «a oC Kaa "Bon a e : 
iid 


A 
2 3 go ce об-е жа T 


& afna tot- ad a^i a) T 


ve alia! a^ 


+a -1) Y 
Е , 


| 
| 
| 
[et ner? кекш “Ky 
| 
| 
E 
| 


= aKes T “Hg 


where 
Ку = Ла - 458. + 22 ваа 
K = 4.48 + Ba 
2 Jom d ё 
2 
к + 450.46, 540: 
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2 
_ A Aa A Ba 
Қз EA E - 
Kk. = та еі) 
5 = 
(а-а +) 
a* a Tes -1) 
6 <% -о 4 n a + 1) ыла “E + 
E. z au 
5. Ж fial 18 „м ү Sia =l) 
SB 2. on? pee Hola) alate (oan) 
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